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Abstract 


We investigate observational constraints on the variable Chaplygin gas model from the 
gold sample of type Ia supernova data and the recent measurements of the X-ray gas mass 
fractions in galaxy clusters. Combining these databases, we obtain a tight constraint on 
the two model parameters. Our results indicate that the original Chaplygin gas model is 
ruled out by the data at 99.7% confidence level. 
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1 Introduction 


Recent observations of type Ia supernovae (SNe Ia) suggest that the expansion of the 
universe is accelerating and that two-thirds of the total energy density exists in a dark 
energy component with negative pressure [I]. In addition, measurements of the cosmic 
microwave background (CMB) and the galaxy power spectrum [3] also indicate the 
existence of the dark energy. The simplest candidate for the dark energy is a cosmological 
constant A, which has pressure Pa = —p,. Specifically, a reliable model should explain why 
the present amount of the dark energy is so small compared with the fundamental scale 
(fine-tuning problem) and why it is comparable with the critical density today (coincidence 
problem). The cosmological constant suffers from these problems. One possible approach 
to construct a viable model for dark energy is to associate it with a slowly evolving and 
spatially homogeneous scalar field ¢, called “quintessence” {4} [5]. Such a model for a broad 
class of potentials can give the energy density converging to its present value for a wide set 
of initial conditions in the past and possess tracker behavior. 


Recently, the Chaplygin gas model was proposed as an alternative to the cosmological 
constant in explaining the accelerating universe [6]. The Chaplygin gas is characterized by 
an exotic equation of state P = —A/p, where A is a positive constant. An attractive feature 
of the model is that it can naturally explain both dark energy and dark matter. The reason 
is that the Chaplygin gas behaves as dust-like matter at early stage and as a cosmological 
constant at later stage. Some possible motivations for this model from the field theory point 
of view are discussed in Refs. [7]. The Chaplygin gas appears as an effective fluid associated 
with d-branes [8] and can also be derived from the Born-Infeld action [9]. An interesting 
range of models was found to be consistent with SN Ia data [10], CMB experiments [I] 
and other observational data [2]. The Chaplygin gas model has been extensively studied 
in the literature [3]. 


However, the Chaplygin gas model produces oscillations or exponential blowup of the 
matter power spectrum that are inconsistent with observation [4]. In Ref. [I5], we con- 
sidered a variable Chaplygin gas (VCG) model and showed that it interpolates between 
a universe dominated by dust and a quiessence-dominated one described by the constant 
equation of state. Furthermore, we showed that the model corresponds to a Born-Infeld 
tachyon action [15]. Recently, the model parameters were constrained using the location of 
peaks of the CMB spectrum and SN Ia data [16]. In this paper we consider observational 
constraints on the VCG model from the gold sample of 157 SN Ia data and the recent mea- 
surements of the X-ray gas mass fractions in 26 galaxy clusters. We perform a combined 
analysis of these databases and obtain the confidence region on the two parameters. Our 
results indicate that the original Chaplygin gas is ruled out at 30 confidence level. 


2 Variable Chaplygin Gas 


Let us now consider the VCG characterized by the equation of state 
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where A(a) is a positive function of the cosmological scale factor a. This assumption is 
reasonable since A(a) is related to the scalar potential if we take the Chaplygin gas as a 
Born-Infeld tachyon field {9} [L5]. In a spatially flat Friedmann-Robertson-Walker (FRW) 

universe, the energy conservation equation is 
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By inserting Eq. (J) into the above equation, one finds that the VCG density evolves as 
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where B is an integration constant. Given a function A(a), Eq. B) allows us to obtain a 
solution p,(a). We assume A(a) is of the form A(a) = Aga~", where Ap and n are constants. 
This ansatz has the following important features: a) the VCG model with n = 0 reduces 
to the original scenario, b) we will see that the VCG behaves as a quiessence rather than a 
cosmological constant at late times and c) the function p(a) can be calculated analytically. 
Then from Eq. (B) it follows that 


=, (4) 


Note that n = 0 corresponds to the original Chaplygin gas scenario, in which the Chaplygin 
gas behaves initially as dust-like matter and later as a cosmological constant [6]. However, 
Eq. (@ shows that, in the VCG scenario, it interpolates between a dust-dominated phase 
and a quiessence-dominated phase described by the constant equation of state w = —1 + 
n/6 [L]. From Eq. (4) we get the present value of the energy density of the VCG 
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where the present value of the scale factor is normalized to unity, i.e., ag = 1. Defining 
B, = B/ Žo, Eq. @ takes the form 


pul2) = po [Bal +2)°+(1-B a+"), (6) 


where z = 1/a — 1 is redshift. In the spatially flat FRW metric the Friedmann equation 


can be written as 
K 


H= z (P + px), (7) 


where H = a/a is the Hubble parameter, K? = 87G is the gravitational coupling and p is 
the energy density of the baronic matter. Substituting Eq. (@) into the Friedmann equation 


(J gives 


H(2)/HE = O4(1+ 2)? + (1-9) 
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where Ho = 100hkms~'Mpc™! is the present value of the Hubble parameter and Q, is the 
density parameter of the baryonic matter component. Then it is straightforward to show 
that the luminosity distance dz, and the angular diameter distance d4 in the spatially flat 
FRW universe are respectively given by 

c z dz 
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dy 


3 Supernova Ia Constraints 


Let us now consider constraints on the VCG model from the gold sample of 157 SN Ia data 
compiled in Ref. [18]. The parameters in the model are determined by minimizing 
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where g; is the total uncertainty in the observation, Hobs is the observed distance modulus 
of SNe Ia, and fmoa(z;) is the theoretical distance modulus 


dr(zi) 
Mpc 


Hmoa (zi) = 5 logig + 25. (12) 
To determine the likelihood of the parameters B, and n, we marginalize the likelihood 
function L = exp(—x?/2) over h and Qs. We adopt Gaussian priors such that h = 0.720.08 
from the Hubble Space Telescope Key Project and Qh? = 0.0214 + 0.0020 from the 
observed abundances of light elements together with primordial nucleosynthesis [20]. The 
results of our analysis for the VCG model are displayed in Fig. I] We show 68.3%, 95.4% 
and 99.7% confidence level contours in the (B,, n) plane. The best-fit model parameters 
and marginalized 1ø error bars are B, = 0.223035 and n = —3.0135 with x2, = 174.253. 
The results show that the two parameter B, and n are highly degenerate. It is interesting to 
note that the dark energy component with w < —1 is favored, which allows the possibility 
that the dark energy density in increasing with time. 


201709.00134v1 


chinaXiv 


Figure 1: Probability contours for B, versus n are shown at lo, 20, and 3a when Qiot = 1. 
These constraints use the gold sample of 157 SN Ia data [I8]. 
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Figure 2: Probability contours for B, versus n are shown at lo, 20, and 30 when Qiot = 1. 


These constraints use the X-ray gas mass fraction in 26 galaxy clusters [25]. 


4 Constraints from the X-ray Gas Mass Fraction of 
Galaxy Clusters 


The matter content of rich clusters of galaxies is thought to provide a fair sample of the 
matter content of the universe as a whole. The observed ratio of the baryonic to total mass 
in clusters should therefore closely match the ratio of the cosmological parameter Qy/Qm, 
where Q, and Qm are the mean baryon and total mass densities of the universe in units of 
the critical density. The combination of robust measurements of the baryonic mass fraction 
in clusters with accurate determinations of Q, from cosmic nucleosynthesis calculations 
can therefore be used to determine Qm BI]. The measurements of the apparent redshift 
dependence of the baryonic mass fraction can also, in principle, be used to constrain the 
geometry and dark energy density of the universe [22]. The first successful application 
of such a test was carried out by Allen et al. using a small sample of X-ray luminous, 
dynamically relaxed clusters with precise mass measurements, spanning the redshift range 
0.1 < z < 0.5 [23] (see also [24]). Recently Allen et al. present a significant extension 
and obtained a tight constraint on the mean matter density and dark energy equation of 
state parameter. The clusters sample is significantly larger and includes 26 X-ray luminous, 
dynamically relaxed systems spanning the redshift range 0.07 < z < 0.9 [25]. We will use 
this database to constraint the VCG model. To determine the confidence region of the 
model parameters, we use the following x? function 
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where o is the measured X-ray gas mass fraction fas with the defaut standard cold 
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dark matter (SCDM) cosmology, Cgas,; is the symmetric root mean square errors, and fgs 


is the model function 
mod 
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where the bias factor b = 0.824 + 0.089 is a parameter motivated by gas dynamical 
simulations, which suggest that the baryon fraction in clusters is slightly depressed with 
respect to the universe as a whole, and the effective matter density parameter is 


Oe = 0, + (1 — %)y/ Bo. (15) 


Adopting Gaussian priors such that h = 0.72 + 0.08 and Q)h? = 0.0214 + 0.0020, the 
68.3%, 95.4% and 99.7% confidence level contours in the (B,, n) plane are shown in Fig. 2] 
The best-fit model parameters and marginalized lo error bars are B, = 0.049+9016 and 
n = 0.5179 with X2im = 24.437. The results favor the original Chaplygin gas model with 


n= 0. 
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Figure 3: Probability contours for B, versus n are shown at lo, 20, and 30 when Qi; = 1. 
These constraints use the gold sample of 157 SN Ia data [IS] and the X-ray gas mass fraction 
in 26 galaxy clusters [25]. 


5 Combined Analysis, Conclusions and Discussions 


Let us now consider a combined analysis of the constraints from the gold sample of SN Ia 
data and the measurements of the X-ray gas mass fraction in galaxy clusters. The fit is 
done by minimizing the following xy? function 


Xai = XSN ag Xray : (16) 


Fig. B] shows the constraints on B, and n obtained from the analysis of the combined the 
SN Ia and the X-ray gas mass fraction data set. We see that SNIa+ fas; data set provides a 
remarkably tight constraint in the (B,, n) plane, with best fit values B, = 0.05510013 (lo 
error bar) and n = 1.70+9:33 (Lo error bar) with X2im = 205.457. We find that the original 
Chaplygin gas model is ruled out by the data at 99.7% confidence level. Fig. Øland Fig. 5] 
show the probability distribution for n and B,, respectively, marginalized over B, and n in 


the spatially flat FRW universe. 


The VCG model, which plays the role of both dark energy and dark matter in terms 
of a single component, deserves to explore its various observational effects [I6]. In this 
paper we have focused our attention on the gold sample of SN Ia data and the recent 
measurements of the X-ray gas mass fractions in galaxy clusters. Adopting simple Gaussian 
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Figure 4: Probability distribution for n marginalized over B, obtained from the analysis of 
the combined SNIa+ fas data set in the VCG model. 
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Figure 5: Probability distribution for B, marginalized over n obtained from the analysis of 
the combined SNIa+ fg; data set in the VCG model. 


priors of h and Q,h?, we have obtained stringent constraints on the two parameters, Bs 
and n, which describe the effective density parameter of the dark matter component and 
the effective equation of state of the dark energy component, respectively. The combined 
analysis of these databases shows that the original Chaplygin gas model is ruled out and the 
universe tends to be quiessence-dominated rather than phantom-dominated at 30 confidence 
level in the VCG scenario. It would be interesting to investigate the evolution of density 
perturbations in this VCG model. 
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